Astronomy & Astrophysics manuscript no. pentericci 


© ESO 2008 


February 5, 2008 





Physical properties of z~ 4 LBGs: differences between galaxies 

with and without Lya emission 



L.Pentericci 1 , A. Grazian 1 , A. Fontana 1 , S. Salimbeni 1 , P. Santini 1 , C. De Santis 1 , S. Gallozzi 1 , and E. Giallongo 



Key words. Galaxies:distances and redshift - Galaxies: evolution - Galaxies: high redshift - Galaxies: fundamental para meters • 



o 
o 

l— i INAF - Osservatorio Astronomico di Roma, Via Frascati 33, 1-00040, Monteporzio, Italy 

S: 

ABSTRACT 

Aims. We have analysed the physical properties of z ~ 4 Lyman Break Galaxies observed in the GOODS-S survey, in order to investigate the 
^ ■ possible differences between galaxies where the Lya is present in emission, and those where the line is absent or in absorption. 

' Methods. The objects have been selected from their optical color and then spectroscopically confirmed by Vanzella et al. (2005). From the 
~- — . public spectra we assessed the nature of the Lya emission and divided the sample into galaxies with Lya in emission and objects without Lya 
f^*) ■ line (i.e. either absent or in absorption). We have then used the complete photometry, from U band to mid infrared from the GOODS-MUSIC 
Cf*i ' database, to study the observational properties of the galaxies, such as UV spectral slopes and optical to mid-infrared colors, and the possible 
, differences between the two samples. 

Finally through standard spectral fitting tecniques we have determined the physical properties of the galaxies, such as total stellar mass, stellar 
' ages and so on, and again we have studied the possible differences between the two samples. 
r~j Results. Our results indicate that LBG with Lya in emission are on average a much younger and less massive population than the LBGs 
O |' without Lya emission. Both populations are forming stars very actively and are relatively dust free, although those with line emission seem to 
' be even less dusty on average. We briefly discuss these results in the context of recent models for the evolution of Lyman break galaxies and 
i, • Lya emitters. 
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H , 1. Introduction LBGs population. Many Lya emitters have now been found 

(e.g. Iye et al. 2006, Ouchi et al. 2004, Fujita et al.2003), and 

In the past years large samples of galaxies have been found up seyeral distant large scale structures or protoc i uste rs have been 

to the highest redshifts (Iye et al. 2006, Kashikawa et al. 2006, discovered (e _ g _ Quchi et aL 20 5 , Venemans et al. 2007). 

Bouwens et al. 2004) using techniques that rely on various Each of the tWQ methods suflfers frQm different selection bias: 

color-selection criteria. Among the various methods, one of the ±e tWQ resulting galaxies population are partia ii y overlapping, 

most efficient is the Lyman break dropout technique (Steidel and k is nQt clear what ig the reMon between them 

& Hamilton 1993) which is sensitive to the presence of the Various scenarios have been proposed t0 exp i ain the prop er- 

912 A break and is effective in finding star forming galaxies. tiefj of Lya emitters Based on rest -f r ame optical photometry 

This method requires a blue spectrum, implying low to moder- of LBGs shapley et al (2001) conc i ude d that LBGs with Lya 

ate dust absorption. It was first designed to select ~ 3 galaxies in emission are w > (ages larger than few xl0 8 yr)> while 

(Steidel et al. 1996, Madau et al. 1996) with a typical magni- » young - (ages less than ^ 100 Myrs) LB Gs have Lya in ab- 

tude limit R>25.5, and then extended to higher redshift (Steidel sorption This could be explained if the young galaxies contain 

et al. 1999; Ouchi et al. 2004; Giavalisco et al. 2004; Dickinson dust which absorbs the Lyff photonSi while the older galaxies 

et al. 2004). are more quiescent with less dust and superwinds which allow 

An alternative technique to find such distant objects is to search tbe EyQ , p botons t0 esca p e 

for Lya emission, through very deep, narrow band imaging In alternativei other groups have sugg ested that strong Lya 

at selected redshift windows, as first shown by Cowie & Hu emiUers m instead young star forming galaxieSi as derived 

(1998). Lya emitters (LAEs) are generally selected to have from the blue colors and high equivalent widths of the Lya 

high restframe Lya equivalent width, typically EW > 20A, emitters (e g Le Fevre et al. 1996, Malhotra & Rhoads 2002; 
with no constraint on the continuum. Therefore, this method & Malhotra 20 01 ; Tapken et al. 2004; Keel et al. 2002). 

tends to select much fainter galaxies, compared to the general Finally some authors haye suggested that galaxies could have 

more than one Lya bright emission phase (e.g. Thommes & 
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Meisenheimer 2005). An initial - primaeval - phase in which 
dust is virtually non-existent, and a later secondary phase in 
which strong galactic winds as observed in some Lyman break 
galaxies facilitate the escape of Ly-a photons after dust has al- 
ready been formed. 

Clearly, it would be important to understand the real relation 
between galaxies with Lya emission and the general LBG pop- 
ulation, so that properties of the overall high redshift galaxy 
population, such as the total stellar mass density, can be better 
constrained. 

To asses this issue we have analised a sample of LBGs selected 
as B-dropouts from the GOODS-South sample, and with VLT 
spectroscopic confirmation (Vanzella et a. 2005, Vanzella et al. 
2006). Given the good quality of the spectra it was possible 
to asses whether the galaxies have Lyo- in emission or the line 
is absent and/or in absorption. We have then analised how the 
observed properties and the derived physical properties such as 
total stellar mass and age depend on the nature of the Lya emis- 
sion. 

All magnitudes are in the AB system (except where otherwise 
stated) and we adopt the A-CDM concordance cosmological 
model (H = 70, Q M = 0.3 and Q A = 0.7). 

2. Sample and observational properties 

From the GOODS-S public data survey we have selected all 
galaxies with spectroscopic redshift in the range 3.4 < z < 4.8. 
These galaxies were initially selected as B-band (or V band for 
those at z > 4.5) dropouts by Giavalisco et al. (2004) (see this 
paper for the color selection criteria adopted) and with a z-band 
constrain z < 26. Spectroscopic observations were carried out 
with FORS2 in the frame of the GOODS project (Vanzella et 
al. 2005, Vanzella et al. 2006). The lower redshift cut is simply 
given by the actual availability of FORS2 spectra, given that the 
Lya line falls at ~ 5300 A, which is at the limit of efficiency for 
the observation-setup used by Vanzella et al. (2006) (FORS2 
plus the 3001 grismo the highest redshift cut was chosen since 
at redshift higher than ~ 4.8 , the spectroscopic confirmations 
are almost exclusively based on the presence of the Lya line in 
emission, while the objects with possible Lya in absorpion (or 
absent) become progressively more difficul to identify. In fact 
the most distant object identified with Lya- in absorption is at 
z=4.788. 

From the public GOODS data we have retrieved the spectra, 
of quality A and B. The spectra with quality flag C were only 
included if our independently measured photometric redshift 
(from 14 bands photometry, Grazian et al. 2006) was in agree- 
ment with the spectroscopic redshift, within the errors. In the 
GOODS survey spectra are classified as emitters or absorbers, 
depending on the nature of the Lya line. After removing AGNs, 
in the above redshift range there are 47 galaxies, of which 19 
have the line in emission, and 28 have Lya in absorption or ab- 
sent. The implied ratio of line emitters/total is quite similar to 
the proportion of Lya emitters that was found at redshift z ~3 
in a large sample of LBGs by Steidel et al. (2000). We will call 

1 At lower redshift spectroscopic observations have been carried out 
with VIMOS but data are still not available to the public. 



the first sample LBGl (LBG with Line emission) and the sec- 
ond sample LBG^ (LBGs with No line). 
In Figure 1 (upper panel) we show the redshift distributions of 
the two samples, which are similar. From our GOODS-MUSIC 
database we have then gathered the complete 14 bands mul- 
ticolor information, extending from the U band to the Spitzer 
8/mi band (for details of the catalog see Grazian et al. 2006). 
We have then studied the observational properties: in particular 
we have derived the restframe luminosity at 1400 A, the UV 
continuum slope B and the z-[3.6] color which enconpasses the 
4000 A break. The average observational properties of the two 
samples, as well as the physical properties that will be derived 
in the next section, are reported in Table 1, together with the 
uncertainties, derived using the standard deviation of the mean. 
In the last column of the Table we report, for each parame- 
ter, the probability value P given by the Kolmogorov-Smirnov 
test (KS-test) which tries to determine if the two datasets differ 
significantly. The KS-test has the advantage of making no as- 
sumption about the distribution of data. 

The U V continuum slope was derived from the i-z color, which 
spans a restframe wavelength range from ~ 1500 A to ~ 1800 A 
(at the mean redshift of the sample z ~ 4, and considering the 
central wavelengths of the ACS filters). Although the restframe 
wavelength baseline is not large (and perhaps a i—J color would 
be preferred), we have used the i - z color for two reasons: for 
a very few galaxies we do not have J band information avail- 
able (see Grazian et al. 2006 for details); second, the UV slope 
derived in such way is more easily comparable to other UV 
slope values in the literature. Like other authors, we assume a 
standard power law spectrum with slope B (/^ oc so that 
a spectrum that is flat in f Y has B = —2). The 1216 A break 
only starts to enter the i filter at z>4.7 so we neglect the impact 
of neutral hydrogen absorption for all galaxies. The measured 
slopes are plotted in Figure 1 (lower panel) and range between 
-1 and -2.5 (with few exeptions). The uncertainties on the de- 
termined values range from 0.2 for the brightest objects to 0.8 
for the fainter ones: in the figure we have indicated a median 
error on the individual values of 0.5 in the upper right corner. 
The average value is (B) = -1.8 + 0.13 for the overall sample, 
and (fi) = -2.0 + 0.11 and -1.7 + 0.13 considering the LBG L 
and LBGn separately. The UV slopes we find are consistent 
with the prediction of models for unobscured continuosly star 
forming galaxies (e.g. Leitherer et al. 1999). They are also con- 
sistent with what found by other authors at a similar redshift: 
for example Overzier et al. (2006), for redshift z=4. 1 LBGs and 
LAEs, find 6 = -1.95 , Venemans et al. (2005) find B = -1.65 
for Lya emitters at z=3.1. The values are also similar to the 
average B = -1 .8 + 0.2 of V606 dropouts in GOODS found by 
Bouwens et al. (2006). 

On the other hand our values are a bit steeper than what found 
by Papovich et al. (2001) for LBGs at z=3 «/?) = - 1 .4). Finally 
Shapley et al. (2003), found less steep values and claim a con- 
siderable difference between the LBGl and LBGn at z~3, with 
the former having a steeper slope, with a positive dependence 
of slope on the Lya equivalent width. 

As can be seen in Figure 1, the values of the UV slopes for 
LBGl and LBGn are largely overlapping and a K-S does not 
reject the hypothesis that the two samples could be drawn from 
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the same underlaying distribution (the probability P is at ~ lcr 
see Table 1). However the median values are indeed somewhat 
different, and the trend is the same as what found by Shapley 
et al. (2003) given that the LBGl are bluer on average than the 
LBGn- The difference between our two samples (AB ~ 0.3) is 
actually similar to that found by Shapley et al. (2003), who have 
a AB = 0.36 between their most extreme groups i.e. the strong 
emitters and the strong absorbers (see Table 3 of the mentioned 
paper), but the significance of their AB is higher, given that they 
have a much larger sample with 800 objects in total, and there- 
fore the uncertainties on the average B values are much lower 
than in our case. In conclusion, although we cannot claim such 
a strong dependence of B on the Lya properties, we do indeed 
find a mild indication that LBGl are bluer than LBGn ■ To 
support this result, a similar difference is also found for LBG 
galaxies by Vanzella et al. (2006), who performed spectroscopy 
of B dropouts from the GOODS sample, and then stacking the 
spectra of galaxies with Lya in emission and in absorption sep- 
arately found that the former have a systematically bluer UV 
continuum than the latter. 

We find no evidence for a slope-magnitude relation as found by 
other authors (e.g. Overzier et al. 2006 ) but the errors on the 
slope are quite large and the wavelength baseline for the slopes 
is rather small. 

In Figure 2 we show the z - [3.6] color plotted versus the to- 
tal [3.6] magnitude: this color encompass the 4000A/Balmer 
breaks, and is sensitive both to dust and stellar ages. The [3.6] 
magnitude corresponds to a rest-frame optical magnitude. The 
strong correlation between the optical flux and color is due to 
the magnitude cut at z = 26 that was used for the sample (see 
dotted line). It is clear from the figure that the LBGn are on 
average both brighter (by ~ 0.7 magnitudes) at restframe op- 
tical wavelengths and redder (by more than half a magnitude) 
than the LBGl ■ Given that all objects are relatively dust free 
(as inferred from the UV spectral slopes and consistently with 
the Lyman Break selection criteria), these differences indicate 
that the LBGn have in general a more evolved stellar popu- 
lation and they are more massive. We will derive the physical 
properties in the next section. 

3. Physical properties 

3.1. Method 

Using the multiwavelength GOODS-MUSIC data and the spec- 
troscopic redshifts, we have then determined the physical prop- 
erties of these galaxies, through a spectral fitting technique. 
The spectral fitting technique adopted here is the same that 
has been developed in previous papers (Fontana et al. 2003, 
Fontana et al. 2006), and similar to those adopted by other 
groups in the literature (e.g. Dickinson et al. 2003, Drory et 
al. 2004). Briefly, it is based on the comparison between the 
observed multicolor distribution of each object and a set of 
templates, computed with standard spectral synthesis models 
(Bruzual & Chariot 2003 in our case), and chosen to broadly 
encompass the variety of star-formation histories, metallici- 
ties and extinction of real galaxies. To compare with previous 
works, we have used the Salpeter IMF, ranging over a set of 
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Fig. 1. Upper panel: the redshift distribution of the two sam- 
ples. Full black line are the LBGl , light dashed line are the 
LBGn ■ Lower panel: the UV slope calculated from the i-z col- 
ors, squares are the LBGl and plus signs are the LBGn ■ 



metallicities (from Z = O.O2Z to Z = 2.5Z ) and dust extinc- 
tion (0 < E(B - V) < 1.1, with a Calzetti or a Small Magellanic 
Cloud extinction curve). Details are given in Table 1 of Fontana 
et al. (2004). For each model of this grid, we have computed 
the expected magnitudes in our filter set, and found the best- 
fitting template with a standard x 1 minimization. The stellar 
mass and other best-fit parameters of the galaxy, like SFR es- 
timated from the UV luminosity and corrected for dust obscu- 
ration (with a typical correction factor of A y ~ 0.4), age, r (the 
star formation e-folding timescale), metallicity and dust extinc- 
tion, are fitted simultaneously to the actual SED of the observed 
galaxy. The derivation of these parameters is explained in de- 
tail in the above paper and in Fontana et al. (2006), where also 
the uncertainties are discussed. In particular we note here that 
the stellar mass generally turns out to be the least sensitive to 
variations in input model assumptions, and the extension of the 
SEDs to the IRAC mid-IR data, tends to reduce considerably 
the formal uncertainties on the derived stellar masses. On the 
other hand, the physical parameter with highest associated un- 
certainty is the metallicity, given that the models are strongly 
degenerate when fitting broad band SEDs; a further limita- 
tion is that only 4 values were allowed in the modelling, i.e. 
0.2,0.02, 1 and 2.5 (respectively subsolar, solar and supersolar 
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Fig. 2. The z-[3.6] color plotted vs. observed flux in [3.6]. 
Squares are the LBGl and plus signs are the LBGn . 



metallicity), but in our samples none of the galaxies turned out 
to have a metallicity of 2.5. 

3.2. Results 

In Figure 3 we plot the distribution of the four basic physical 
parameters for the two samples, namely the total stellar mass, 
the star formation rate, the derived stellar age and the E(B-V). 
In each figure we also indicate with an arrow the average values 
calculated for each sample. 

Both the total stellar masses and the median ages are 
considerably lower for the LBGl compared to the LBGn ■ 
The average mass is (5.0 + 1) x 1O 9 M for the LBG L and 
(2.3 + 0.8) x 1O 1O M for the others, i.e. a factor of almost five 
higher. The K-S test gives a very low value, implying that the 
two populations are different from each other with > 99.8% 
probability. We can therefore conclude that the LBGl are a 
factor of almost 5 less massive than the LBGn ■ 
The median ages are also quite different, with an average of 
200 + 50 Myr for the LBGl , an age distribution that is very 
peaked towards small ages and is basically confined to values 
below 300 Myr. The LBGn on the other hand have an average 
age of 410 + 70 Myr (a factor >2 higher) and there are several 
galaxies with ages exceeding 1 Gyr, which is a considerable 
fraction of the cosmic time at redshift ~ 4. Again performing 
a K-S test, the two populations are different with a probabil- 
ity > 98% . We can conclude that the LBGl are significantly 
younger galaxies compared to the LBGn 
The star formation rates do not differ much for the two sam- 
ples, with median values that are almost equal, (SFR) = 76 + 
34M Q yr~ l and 74 + 22M Q yr~ l respectively and very similar 
distributions. The exctinction values are slightly higher for the 
LBG N with (E(B - V)> = 0.10 + 0.02, while the LBG L have 
(E(B - V)) = 0.07 ± 0.01. In both cases the K-S test does not 
reject the hypothesis that the 2 groups are drawn from the same 
population of galaxies, although the average E(B-V) is higher 
for the LBGn ■ Therefore both LBGl and LBGn are galaxies 



with relatively little dust content, and significant star formation 
rates. In particular the average E(B-V) is quite lower than the 
average values found by Shapley et al (2003) at redshift ~3 and 
Ouchi et al. (2004) at redshift ~4 but we find the same tendency 
for LBGn to have somewhat higher dust content compared to 
LBGl which is also supported by the difference in the UV 
spectral slopes discussed in the previous section. 
The final two physical parameters that come as output from 
the modelling are the metallicity, and r i.e. the star formation 
e-folding timescale. Each of them suffers from larger uncer- 
tainties as compared to e.g. the total stellar mass, as already 
discussed before. The mean values and the total distribution 
for these two parameters are very similar both for LBGn and 
LBGl, with the K-S test indicating that they are drawn from 
the same population. A more interesting and better constrained 
quantity to determine is the age It parameter that is an indica- 
tion of the evolutionary state of the galaxies. It can be shown 
that passively evolving galaxies can be selected according to 
the physical criterion age/r > 4 (Grazian et al. 2006). All 
galaxies in our sample, both LBGn and LBGl have age/r < 4 
as expected: however the difference in age discussed before is 
reflected also on the difference in age/r, which is relatively 
larger for LBGn with a mean value of {age / t^lbgn = 0-99 
while for LBGl the mean value is 0.33. Furthermore all but 
one LBGl have age/r < 1 while amongst LBGn there are 7 
objects with values between 1 and 4. 

From this data we conclude that the LBGl seem to be a 
population of younger and less massive galaxies. Given that 
all galaxies have relatively low dust content, there are some 
indication that LBGl might be relatively less dusty than LBGn 
in general, based on the results presented here and by other 
authors (Vanzella et al. 2006, Shapley et al. 2003 and Overzier 
et al. 2006). The star formation rates are high in all cases and 
are independent on the nature of the Lya emission. 

4. Discussion and Conclusions 

The results obtained here seem to point towards a scenario 
where LBGl are young and relatively small galaxies compared 
to the general LBGs population, are basically unobscured and 
are forming stars at very high rates. LBGn on the other hand, 
are more massive, with masses as high as 1O U M , and span 
a larger range of ages from few tens of Myrs to more than 1 
Gyr. They are forming stars at comparable rates and are also 
relatively unobscured galaxies, although they appear slightly 
redder than LBGl ■ This is therefore partially at variance with 
the only comparable large study, carried out at redshift ~3 by 
Shapley et al. (2003), who found a tendency for LBGl to have 
in general higher ages, lower values of E(B-V) and higher star 
formation rates. A very similar work by Iwata et al. (2005) on 
z~ 3 LBGs from the Hubble Deep Field-South, also examined 
the relationship between spectroscopic features and parameters 
of the best-fit stellar population models in a sample of spectro- 
scopically confirmed galaxies. They did not find any clear cor- 
relation between the Lyo- emission, Lya equivalent widths and 
the results of SED fitting. However the absence of any clear 
trend may be attributed to the very small number of sample 
galaxies in their study (13 in total). 
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Fig. 3. The physical properties of B and V-dropouts. Top left: distribution of total stellar mass; top right the exctinction parameter 
E(B-V) derived from the spectral fitting; middle left the total star formation rate derived from the spectral fitting; middle right: 
stellar ages; bottom left: t the star formation e-folding time-scale; bottom right: the metallicity. Solid black and dashed light 
histograms denote LBGl and LBGn respectively. The blck and light arrows indicate the mean values of each sample. 



Some other observational studies seem to support our results: 
for example Overzier et al. (2006) find similar trends for a small 
sample of LBGs and LAEe in a z=4.1 protocluster, with indi- 
cations that LAEs are young, dust free objects with blue UV 
optical colors and possibly less massive that UV selected LBGs 
with comparable star formation rates. 

At higher redshift, Lai et al. (2006) constrained the ages of Ly<? 
emitters at redshift 5 and found typical ages 50-100 Myr, i.e. 
comparable to our LBGl ■ Finally at redshift 3, Gawiser et al. 
(2006) analised a sample of narrow band selected Lyor emitters 
at z~ 3 and concluded that they appear to have much less dust 
and lower stellar masses, compared to the generic LBG popula- 
tion at the same redshift. They were not able to constrain ages, 
which had a large range of allowed values. However given that 
they started from a narrow band sample, their objects are in 
general much fainter than the general LBGs population at the 
same redshift and clearly they are not selected to have the same 
broadband optical colors as the LBGs at the same redshift. In 
our case the optical colors and absolute restframe UV luminos- 
ity of the LBGn and LBGl are equal. 

Our results give support to models that consider the Lya 
emitters as precursors of spheroid galaxies, observed during 
the initial starburst phase, which is confined to a short period 
after its onset, due to the rapid formation of dust. Interestingly 
one such model (Thommes & Meisenheimer 2005) attempts 
to constrain the duration of this initial phase from observa- 



tional data, finding T^ya ^ 350 Myr. A very similar results 
was recently found by Mori & Umemura (2006) from simula- 
tion with a N body/hydrodynamics code: they follow the early 
stages of galaxy formation and find a first phase characterised 
by strong Lya emission, which after 3x10 s years quickly de- 
clines to fluxes below the observable level. The ages found in 
these models, are very similar to the range of allowed ages for 
our LBGl as shown in Figure 3; on the other hand many of 
the LBGn have ages exceeding 300 Myr. Clearly there are 
also very few LBGl with larger ages and most of all, there 
are several LBGn which are very young, with ages of a few 
tens of Myrs. The first few outliers could be easily accounted 
for given that the modelled physical parameter are obviously 
subject to uncertainties. However the very young LBGn i e. 
primeval objects without Lya emission are hard to fit within 
the simple scenario proposed above. Probably more complex 
models in which the time scales for dust formation in primeval 
galaxies is variable and depends on other factors have to be in- 
voked. In this context, several scenarios have been proposed, 
for example where the differential attenuation of the Lya and 
continuum photons depends on the clumpyness of the medium 
(e.g. Hansen & Oh 2006), or where the timescales for the Lyor 
emission varies for galaxies residing in halos of different sizes, 
in the sense that in the most massive halos the Lyo- luminosity 
declines abruptly after a much shorter time (10 7 ys) than in less 
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Table 1. Average properties of LBGi and LBG^ 



Property 


LBG L 




P(K-S) 




19 


28 




z 


24.89±0.13 


24.85+0.10 


0.99 


[3.6] 


24.6±0.2 


23.8+0.15 


0.02 


P 


-2.0±0.11 


-1.7+0.12 


0.32 


L(1400) (cgs) 


l.lxlO 29 


l.OxlO 29 


0.65 


SFR(M yf--') 


76±34 


74+22 


0.90 


Age (Myrs) 


200±50 


410+70 


0.02 


E(B-V) 


0.07±0.01 


0.10+0.02 


0.58 


Mass (M ) 


5 ± 1 x Iff 


2.3 + 0.8 x 10 10 


<0.002 


Metallicity 


0.55+0.10 


0.48+0.08 


1.0 


T 


4.5+1.4 


5.5+1.3 


0.96 


Age/r 


0.33+0.13 


0.99+0.23 


0.06 



Average observational and physical properties for the samples of 
LBG L and LBG N separately: B is the UV slope derived from the i-z 
color, L(1400) is the total luminosity at 1400Aresframe; SFR is the 
total instantaneous star formation rate; Mass in the total stellar mass 
assembled in the galaxies; t is the star formation e-folding timescale. 
" We remind that this parameter is the most uncertain from the SED 
fit output, and that only four values of metallicity were considered in 
the fitting procedure, see text for more details. 

massive halos (Mao et al. 2007). However a full discussion of 
the models is beyond the scope of this paper. 

Finally it also remains unclear whether the LBGn are ac- 
tually dustier than the LBGl , which would be a natural conse- 
quence of this model: we do indeed find that LBGl are slightly 
bluer than LBG^ (with the slope B determined directly from 
the observed colors) and consequently have on average lower 
E(B-V) parameter, as derived from the spectral fitting. However 
the differences are not very large. An important independent 
check would be to assess the dust content e.g. from far IR ob- 
servations, rather than through the E(B-V) which is subject to 
larger uncertainties and is strongly degenerate with the fitted 
age. 
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